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Abstract 

It is demonstrated that there is a dynamic isospin breaking effect 
in the near threshold 7*iV — > ttN reaction due to the mass difference 
of the up and down quarks, which also causes isospin breaking in the 
ttN system. The photopion reaction is affected through final state 
irN interactions (formally implemented by unitarity and time reversal 
invariance) . It is also demonstrated that the near threshold 7iV — > ttN 
reaction is a practical reaction to measure isospin breaking in the ttN 
system, which was first predicted by Weinberg about 20 years ago but 
has never been experimentally tested. 

Since the discovery of a large mass difference between the up and down 
quarks (m u ~ 5MeV, m d ~ 9MeV) there has been considerable interest in 
the possibility of observing dynamical isospin breaking in the pion-nucleon 
system Jj], |2], J|, [J. The theoretical consensus is that the magnitude of 
isospin breaking is not (m d — m u ) / {m d + m u ) ~ 27%, but instead is ~ {m d — 
m u )/A qcd ~ 2% 0, where A qcd ~ 0.2 GeV. In this paper it will be shown for 
the first time that the 7^ — ► ir°p reaction is an excellent candidate to measure 
dynamic isospin violations due to the up, down quark mass difference. 

Weinberg first showed that there is an isospin violating effect in the s wave 
71N scattering length a(irN) due to the up, down quark mass difference 
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This predicted effect, which occurs for tt°N scattering or charge exchange 
but not for n^N scattering, has more recently been calculated by Chiral 
Perturbation Theory (ChPT) [|], The predicted magnitude of this effect 
is the same (to within a factor of a/2) in ir°N scattering and charge exchange 
reactions. However, since the magnitude of a(Tr°N) is small, the relative 
magnitude of the isospin violating term is ~ 30 % ^ f|. By contrast, for 
charge exchange where the isospin conserving amplitude is larger, the relative 
isospin violation is estimated to be ~ 2 to 3 % ]1], f|. 

In this paper the (dynamic) isospin breaking effect of the up and down 
quark mass difference is shown to be present in the near threshold 7 N — ■> 7rN 
reaction. This is in addition to the well known (static) isospin breaking 
effect due to the threshold difference between the jp — > n°p and 7p — > 7i + n 
channels. Of this ~ 2/3 is due to the Coulomb effect and ~ 1/3 is caused 
by the up, down quark mass difference. The separation of the two threshold 
energies leads to the prediction of a unitary cusp in the 7p — > rr°p reaction 
near the u + n threshold due to the two step ^p — > n + n — > 7T°p reaction 
II 0> H> HI- The magnitude of this unitary cusp is due to the large ratio of 
the electric dipole multipole amplitudes, R = E Q+ (jp — > n + n)/ E 0+ (jp — > 
7i°p) ~ —20 0. The first derivations of the unitary cusp used a K matrix 
approach to calculate the effect of charge exchange in the final state. Due to 
their single scattering approximation, these calculations violate unitarity as 
will be discussed below. The chiral perturbation theory (ChPT) calculations 
are basically isospin conserving but the biggest isospin non-conserving effect 
due to the pion mass difference (which is almost purely electromagnetic in 
origin) is inserted by hand ||. These calculations also violate unitarity due 
to their truncation at the one loop (single rescattering) level. There has been 
recent progress in performing ChPT calculations in nN scattering which take 
the isospin violations due to the Coulomb interaction into account ||. 

The approximations made in the K matrix [0] and published ChPT pa- 
pers can be overcome by using a 3 channel S matrix approach in which 
unitarity and time reversal invariance are satisfied (a preliminary version of 
this work has been presented previously 0). This has the advantage that 
both static and dynamic isospin breaking, and final state multiple scattering 
(to all orders), are taken into account. The S matrix for the 3 open channels 
(7P, n°p, n + n) can be written as : 
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where S 7 = e 2l5 ~> , S = cos (pe 2l5 ° , S c = coscj)e 2l5c represent elastic jp, 7i°p, 
n + n scattering with phase shifts <5 7 , S , and 5 C respectively, S 0c = S c0 = 
isin^e^' 5o+5c ' ) represents the n + n «-> ir°p charge exchange amplitude, is a 
real number, and cos0 represents the inelasticity due to charge exchange. 
For convenience, the off diagonal matrix elements for the photo-production 
of the ii p and n + n channels are written as iM' Q and iM' c , where Mq and 
M' c are proportional to the multipole amplitudes. For the important case of 
production of s wave pions these are the E +(L +) transverse (longitudinal) 
multipoles. Although not explicitly written here, the S matrix elements are 
for a fixed value of W, the total CM energy, and the quantum numbers 1 and 
j, the ttN orbital and total angular momenta. The notation of using (c) for 
the neutral (charged) channel is conveniently generalized to neutron targets 
where the three channels are jn, n°n, and n~p. 

Time reversal invariance requires that the S matrix be symmetric (Sij = 
Sji) and unitarity requires that S + S = SS + = 1. The form of the 2x2 tcN 
portion of the S matrix has been chosen to be separately unitary and time 
reversal invariant. Eq. |TJ is symmetric by construction. Applying the unitary 
constraint and assuming the weakness of the electromagnetic interaction by 
dropping terms of order e 2 , one obtains [|10|] : 



where A' and A' c are smoothly varying, real functions, of the CM energy, and 
can be identified as proportional to the multipole matrix elements Mo (for 
7P — > ir°p) and M c (for ^p — ► 7r + n) in the absence of the final state ttN 
interaction. These equations show the connection between electromagnetic 
pion production and nN interactions. Note that the phase shifts <5 7 , So, and 
S c appear as e tS in photopion reactions compared to e 2tS in elastic scattering 
because the nN^N) interactions take place in the final (initial) state only. 
Eqs. |I] and |2| are valid for both photo- and electro-production, i.e. for both 
real and virtual photons. In general Mq, M' c , A' q , and A' c are functions of q 2 
(the invariant four momentum transfer) as well as W, 1, and j. To order e 2 
the tiN sector parameters 5q, 8 c , and are functions of W, 1, and j only. 
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In this treatment terms of order e 2 have been neglected with the exception 
of S y which is small compared to both S and S c (except at the 7r° threshold) , 
and the mass difference between the charged and neutral pions which is put 
in by hand. This is the largest 0(e 2 ) effect found in tcN scattering || and 
in photoproduction ||. The effect of the 0(e 2 ) terms have been worked 
out, although not presented here. In general they produce a small coupling 
between the electromagnetic and ttN sectors, so that e.g. the 2x2 ttN sector 
of the S matrix is no longer separately unitary. When the experimental 
situation reaches a sufficient level of accuracy these corrections should then 
be implemented. 

Eq. pi is a generalization of the final state interaction theorem of Fermi 



and Watson ]TTj] who first pointed out the connection between photo-pion 
production and ttN scattering. In their derivation, as in this one, time re- 
versal and unitarity (to order e 2 ) were assumed. However they made the 
additional assumption that if one could neglect e 2 (and higher order terms) 
isospin would be conserved (this was before the time of quarks). This re- 
duces the dimensionality of the S matrix to 2x2 which gives the simple 
solution S(jp — ► TtN : I) — ±\S\e i5 ( nNd) where 1=1/2,3/2 is the isospin 
of the ttN system. If isospin were conserved, so that the two thresholds 
are degenerate, Eq. ||] reduces to the Fermi- Watson theorem. For the near 
threshold region in which the CM pion momentum q — > 0, the s wave phase 
shifts are 5 = [2a 3 + ai]g/3 = (V^g, 5 C = [a 3 + 2ai]g/3 = a^+nq, and 
4>/2 = [ai — a 3 ]q/3 = a cex q where cij is the s wave ttN scattering lengths 
in the isospin 21 = 1, 3, or designated charge states. However isospin conser- 
vation is badly violated in the threshold region due to the threshold energy 
difference between the n + n and 7c°p channels and the additional dynamic 
isospin violating effect due to the up and down quark mass difference. This 
generalization of the Fermi- Watson theorem removes the approximation of 
isospin conservation. Note the interesting feature that below the ir + n thresh- 
old there are only two open channels and Eq. ^ reduces to the form of the 
Fermi- Watson theorem, i.e. M' Q = real number • e lS ° or, equivalently, that 
tan(5 ) = lmE Q+ /ReE 0+ . 

It is of interest to show the connection between the S matrix formula- 
tion which takes the final state scattering into account to all orders and the 
original K matrix derivations which only took a single final state scat- 
tering into account. This single scattering limit can be taken by expanding 
e "5 ~ ]_ _|_ i§ 0j an d neglecting second order terms like Sq5 c , thus obtain- 
ing M ~ A + iA 5 + iA c a cex q c . By examing the region below the n + n 
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threshold where q c is continued as i\q c \, it is observed that tan(5o) ~ Sq ^ 
IhlE'o+ /Re-E'o+ This means that the single scattering approximation does not 
obey unitarity, although in this particular case the numerical effects are not 
large. The ChPT calculations to one loop have a similar unitarity problem 
as was discussed by the ChPT authors ||. 

The connection between ML, to the multipole amplitudes M 0(c ) is ML c s = 
2 v /g 7 g ( c )M ( c ). Note that the multipoles (M) have the dimensions of length, 
while the S matrix elements (M') are dimensionless. The s wave photopro- 
duction cross section is 0o( c )(0) = \MLA 2 /Aq 2 = (qo( c ) / q-y)\Mo^\ 2 . Similarly 
the s wave n + n — > 7r°p charge exchange reaction cross section is cr cex (6) = 
sin 2 0/4g2 = (q /q c )\f cex \ 2 from which one obtains sin0 = 2 y /q^\f cex \ jl2| 



Defining A' Q ^ = 2 v /g 7 g (c)^4o(c) one can then write Eq. | in terms of multipole 
amplitudes as: 

M =e i ^ +<5 °)[A)Cosf + tA c q c \f cex \/cosl] 

Mc =^^+^^008 I + iA q \f cex \/ cos I) U 

In the threshold region Eqs. |2] and |3] lead to a unitary cusp. In this 
case, -> a cea: g c ,cos(f ) — > 1, \ f cex \ — >■ a cex , and A c ~ E 0+ (-yp -> 7r + n) for 
the transverse multipole, or L +(7P — > 7r + n)) for the longitudinal multipole. 
Below the n + n threshold q c is continued as i\q c \. Thus the contribution of 
the two step jp — > it + n — > n°p reaction is in ReM (lmE 0+ ) below (above) 
the 7r + n threshold. The magnitude of the unitary cusp, (3, is defined as the 
coefficient of q c : 



(3 = A c f cex {Ti + n -> 7r°p) /cos(4>/2) ~ -E +(7P -> vr + n)a ce:c (7r + n -> vr>) 

(4) 

le value of j3 can be calculated from the experimental value of a cex {j[ p — > 
Assuming isospin is conserved a(7f + n <-> 7T°p) = —a(-K~p <-» 7r°n) 
IH. For £ 0+ ( 7 p -> Tr+n) the ChPT prediction of 28.2 ±0.6@| is used. 
This is in agreement with the preliminary value of 27.6 ± 0.3 from a recent 
Combining the two experimental (theoretical) numbers 
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experiment 
gives (3 = 3.59 ± 0.17 (J3 = 3.51 ± 0.22) |TJj . 

Although it has not been previously noticed, the same physical mechanism 
that causes the isospin breaking in ttN scattering and charge exchange also 
enters into f3. This follows from unitarity and time reversal (Eqs. [| and [3D, 
combined with the prediction of isospin breaking for the s wave scattering 



5 



lengths for charge exchange reactions Sa cex , that there is an isospin violating 
contribution to ft, 

5ft = E 0+ (jp -> n + n)5a cex (-K + n -> ir°p) , , 

= 5a cex (7T + n -> Ti°p) / a cex (ix + n -> 7r°p) ' ' 

In the second line it has been assumed that there is no quark mass effect on 
E 0+ (jp — ► n + n). It is estimated that Sa cex /a cex ~ 2 to 3% 0, |J. A better 
theoretical calculation of this effect should be performed. 

The electric dipole amplitude Eo + (jp — > 7r°p) with its unitary cusp is 
presented in Figs. [I] and 0. The recent Mainz/TAPS (T7j and Saskatoon 



I8| results for KeE 0+ are presented in Fig. [I], where only the statistical 
errors are shown. The small deviations between the data sets suggest the 
magnitude of the systematic errors. Considering this, the agreement between 
the calculated curves and the data is satisfactory. The curves are the ChPT 
calculation with three empirical low energy parameters used in fitting the 



data ||, and a unitary fit to the data [|17[] using Eq. y which has ft = 3.67 fl4 
and a linear function of photon energy for A with two parameters which were 
fit to the data. For this case R = E 0+ (jp — > n + n)/E 0+ ('yp — > tt°p) ~ —20 
H and the effect of the two step charge exchange reaction is dramatic in 
the real (imaginary) part below(above) the 7r + n threshold. Approximately 1 
MeV above the 7r + n threshold the |Im£?o+| = |Rei?o+|- The negative sign in 
R makes IhlE + have the opposite sign from KeE 0+ . 

It is interesting to note that both the power counting rules of ChPT 
and the constraints of unitarity lead to pion rescattering in the final state 
as a critical dynamical ingredient. Although the unitary and ChPT curves 
both agree with the data for Rei?o+ there is an important difference between 
them. The value of ft = 2.78 || [l4j calculated for the 7p — > ir p reaction 
is smaller then the one calculated using the separately predicted values of 
Eq + {^p — > 7r + n) and a cex {ji + n — > 7i°p) of 3.51 ± 0.22. This difference can be 
clearly seen in Fig. |2|. The reason for this discrepancy, which was discussed by 
the ChPT authors 0, is due to the fact that the ChPT calculation is carried 
out to one loop which is not sufficient for Im.E'o + . This is a general feature 
of chiral perturbation theory in which the imaginary part of the amplitude 
is not calculated as accurately as the real part, and thus unitarity is only 



approximately satisfied at a given order | 19| . As will be discussed below 
the difference between the unitary and one loop ChPT values of ft can be 
observed in future experiments in which lmE 0+ will be measured directly. 
To accurately measure the magnitude of the unitary cusp and to exploit 
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Figure 1: KeE 0+ (in units of 10~ 3 /m 7r ) for the jp — > ir°p reaction versus 
photon energy k. The dashed dot curve is the ChPT fit || and the solid 
curve is the unitary fit (Eq. [T]) JT7J. The solid diamonds (circles) are the 
Mainz (Saskatoon) points |T7| , |I8| . The errors are statistical only. 
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Figure 2: lmE 0+ (in units of 10~ 3 /m 7r ) for the 7p — > ir°p reaction versus 
photon energy k. The dashed dot curve is the ChPT calculation || and the 
solid curve is the unitary calculation (Eq. |T]). 
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the connection between electromagnetic pion production and low energy irN 
interactions, one must measure lmE 0+ . In photoproduction this requires 
experiments with polarized beams and/or targets ||. For the sake of brevity 
the formulas connecting the cross sections and polarization observables to the 
multipoles [^(| will not be quoted here. To briefly demonstrate the power of 
polarized photo-pion experiments, two asymmetries are shown in Fig. ^: £ for 
linearly polarized photons with an unpolarized target; and T for unpolarized 
photons but with a target polarized normal to the reaction plane. The results 
presented in Fig. || use the p wave predictions of chiral perturbation theory 
and the unitary fit to E Q+ discussed above |T7[]. S is primarily sensitive 
to the p wave multipoles and since the unitary fit has essentially the same 
p wave multipoles as ChPT, the curves for ChPT and the unitary fit are 
almost identical. By contrast, T is sensitive to a linear combination of p wave 
multipoles times IhlEo+, and shows its rapid rise above the n + n threshold. 
For T, the large difference in IulEo+ between the unitary fit |17| and ChPT 



to one loop || should be straightforward to distinguish experimentally. A 
proposed experiment with tagged photons using an active, polarized proton 
target at Mainz estimates that (3 can be measured to ~ 1 to 2% Similar 
results could be obtained using a laser backscattering source |22|, ^3|] or small 
angle electron scattering with polarized, internal, targets in a storage ring 
facility |M . 



There are several possible strategies to test isospin conservation from 
measurements in the near threshold 'jN — > tcN reaction. A definitive demon- 
stration of this effect would involve precision measurements of the real and 
imaginary parts of Eq + for the unitary cusp contributions to the jn — > 7i°n 
and 7p — > n°p reactions. Although the predicted effect of ~ 2 to 3% in (3 
is larger than the usual ~1% effect expected on the basis of electromagnetic 
effects, it represents a serious experimental challenge to accomplish this goal, 
particularly since free neutron targets are not available to make this measure- 
ment. A simpler strategy would be to perform a precision measurement of the 
unitary cusp for the 7p — > n°p reaction, extract the value of (3, and compare 
the result of with the value obtained using unitarity and isospin conservation 
quoted above. An equivalent strategy is to obtain a cex {ix + n — > ir°p) from 
the measured value of (3 and compare it to the accurately measured value 
of a cex (ir~p — > 7T°n) [[Tj|, obtained from the width of the Is state in pionic 
hydrogen, to see if there is any deviation from the pure isospin prediction 
that these values are equal and opposite It is important to note that 
there is a PSI proposal |13] to improve the pionic hydrogen measurement so 
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Figure 3: The polarized photon (£) and polarized target (T) asymmetries 
(in %) at a tt° center of mass angle of 90° for the jp — > ir°p reaction versus 
photon energy k. The lower curve is the ChPT calculation |J for E. The 
two curves for T which rise rapidly at the n + n threshold at 151.4 MeV are 
the ChPT calculation (dashed curve) and the unitary calculation (solid 
curve). 
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that the error in a cex {ii~p — > n°n) will be reduced to ~ 1%. 

It is therefore of interest to discuss the information that might be ob- 
tained from measurements of the s wave ttN scattering lengths, a(nN), from 
photopion reactions with the values that have been obtained using conven- 
tional pion beams. The values for a(7rN) for several physical channels are 
presented in Table [I] (the appropriate isospin relations |^5[ were used when 
required). The most accurate and direct measurements come from pionic hy- 



drogen experiments JT3[. In order to test isospin symmetry the PSI group has 
also measured pionic deuterium in order to get at the n~n scattering length 
]l3fl . Unfortunately the two body corrections in deuterium are large and add 
significantly to the error which makes the isospin conservation test uncertain 
at the required level of accuracy, so that these values are not quoted here. 
The results of the two most recent empirical analyses of the irN data for 
pion kinetic energies greater than ~ 30 MeV, extrapolated to threshold, are 
also presented |27|, |26| in Table [I]. In the SAID study the scattering lengths 
come from a dispersive analysis, the errors are hard to ascertain and are not 
quoted [p^| . The chiral perturbation theory results || are consistent with 
experiment. If isospin symmetry is exact there are two independent scatter- 
ing lengths in the I =1/2, 3/2 states, and six possible physical 7rN elastic 
scattering and charge exchange reactions, so that there are many possible 
tests if three or more channels can be measured precisely. This is precisely 
the gap that can be filled by measurement of the final state 7rN interactions 
in photopion production, e.g. ir°p or n + n — > 7T°p in the 7p — > n p reaction. 
Another test of isospin conservation has been made for medium energy 



7rN scattering (pion kinetic energy from 30 to 100 MeV) |27| , |2q| . In both 
cases the scattering amplitudes in the 1=1/2 and 3/2 states were obtained 
from the data for ir ± p elastic scattering. From these, the prediction for the 
n~p — > 7T°n charge exchange amplitude made on the assumption of isospin 
conservation was compared to the empirical charge exchange amplitude. An 
isospin violation at the ~ 7% level has been found by both analyses [E7|, |28[ , 
primarily in the s wave amplitude. These analyses depend on the quality of 
the data and on the accuracy of the Coulomb corrections, which have been 
criticized as being inconsistent with the strong interaction calculations that 
were employed || . If isospin is indeed violated at the 7% level, it is not clear 
how to relate this to the isospin breaking predictions for the s wave scattering 
length |], f|] which is applicable at lower energies. What is required is an 
extension of the isospin breaking calculations to medium energies. A test of 
the existing predictions requires experiments at lower energies which would 
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Table 1: S-wave ttN scattering lengths a(irN) for several channels in units 
of 10~ 2 /m 7r (N = n or p) 



Channel 


Pionic 
Atoms 
11 


Matsinos 

m 


SAID 

(Sp98) 

in 


ChPT 
§ 


Estimated 
isospin 
breaking ||, f| 


TT~p 


8.83(0.08) 


8.14(0.10) 


8.83 


8.70(0.86) 


~ 


7T~P — > 7T°n 


-13.01(0.59) 


-10.93(0.08) 


-12.5 


-12.45(0.75) 


~ 2 to 3 % 


TT°N 




0.41(0.09) 


0.0 


-0.1(0.7) 


~ 30% 



measure the s wave scattering lengths. 

In addition to the measurement of lmE 0+ above the n + n threshold dis- 
cussed above, one could contemplate the more difficult measurement when 
only the ir p channel is open. As was discussed previously, from a measure- 
ment of both the Re and Im parts of E 0+ (or L 0+ ) one obtains 5q, for which 
there is a predicted ~ 30% isospin breaking effect. However, because of the 
small expected size of 5q ~ 1° in the energy region below the ir + n threshold, 
this is a very difficult task. It may be at the limits of feasibility using either a 
laser backscattering source, where an intense photon flux is concentrated in a 
small energy interval [^, |23|| , or with internal target, small angle scattering 
with polarized internal targets 34" . 

In conclusion, the connection between the j*p — > 7r°p reaction and ttN 
scattering has been presented in a rigorous, model-independent way, which 
is unitary and time reversal invariant, and where the isospin breaking due 
to the threshold difference between the n°p and ir + n channels, and the mass 
difference between the up and down quarks, is taken into account. This 
leads to a predicted unitary cusp due to the two step •yp — > n + n — > ir°p 
charge exchange reaction. The magnitude of the unitary cusp is given by 
(3 = E 0+ ('jp — > n + n) ■ a cex (ir + n — > 7r°p). This unitary cusp has been recently 
observed in photoproduction experiments |T8|] . It has been shown for the 
first time that there is a dynamical isospin breaking effect in the value of /3, 
due to the mass difference of the up and down quarks, in electromagnetic pion 
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production. This is linked by unitarity and time reversal invariance to a pre- 
dicted quark mass effect in n°N scattering and pion charge exchange [5], || ■ 
At present there are accurate measurements for a n - p and a cex {ji ~p — > ir°n) 
from pionic atoms performed at PSI [13|]. To check isospin conservation re- 
quires at least one more precision measurement in another charge channel. 
This is more readily performed in electromagnetic meson production for two 
reasons: 1) in order to accurately measure the s wave scattering length, it is 
important to work at very low energies (e.g. < 10 MeV) at which irN experi- 
ments are hard to perform since the low energy charged pions decay and also 
since 7r° beams can not be made at any energy; and 2) in electromagnetic 
pion production one can access charge states that cannot be reached with 
conventional pion beams (this is important for isospin checks). It is shown 
that photoproduction experiments with polarized targets can lead to a pre- 
cise measurement of the small but interesting isospin violating effects. In 
addition to the polarization observables in photoproduction discussed above, 
similar information can be obtained in threshold 7r° electroproduction from 
the combination of the transverse- longitudinal (TL and TL') structure func- 



tions |^0|. This possibility will be discussed in a future publication. Finally 
we stress that the important program of precisely measuring the irN and 
photopion reactions that are required to test the predicted isospin violation, 
is difficult but feasible. 
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